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We report real space imaging measurements of inelastic Friedel oscillations. The inelastic elec- 
tron tunneling spectroscopy, using scanning tunneling microscopy, around dimers of dichlorobenze 
adsorbates on Au(lll) surface display clear spatial modulations that we attribute to inelastic 
scattering at the molecular sites caused by molecular vibrations. Due to local interactions between 
the adsorbate and the surface states, the molecular vibrations generate a redistribution of the 
charge density at energies in a narrow range around the inelastic mode. Our experimental findings 
are supported by theoretical arguments. 
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Vibrational modes are ubiquitous to molecular struc- 
tures in disciplines ranging from molecular electronics 
and chemical analysis to biology. Inelastic scattering at 
these vibrational modes plays an essential role for the 
properties of most materials, including their electric and 
thermal conductivities, e.g. of metallic clusters [TJ [2], 
molecules [2H5|, superconductors [6], surfaces [7J[8|, car- 
bon nanotubes and graphene [9j[10]. Experimentally, in- 
elastic transitions are probed by measuring the energy 
loss required to trigger the inelastic mode. Here, we 
present for the first time real space images of inelastic 
standing waves around adsorbed molecules. These waves 
result from a superposition of surface state electrons be- 
fore and after scattering inelastically at the adsorbate. 
In analogy to charge density waves arising from elastic 
scattering at defects (Friedel-like oscillations) [TTJ [12], 
we name these waves inelastic Friedel-like oscillations. 
Apart from the inherent beauty of such patterns, we an- 
ticipate that investigations of inelastic Friedel oscillations 
will provide a better understanding of vibrational exci- 
tations, coupling of electrons to molecules, and their in- 
fluence on material properties. 

In recent years the ability to probe quantum mechan- 
ical properties of e.g. single atoms [TJ [2] and molecules 
[3l [4| using scanning tunneling microscopy (STM) has 
become remarkable. The mesurements are stretching 
towards atomic resolution even of collective modes as 
phonon excitations in metallic surfaces [7 and their 
transmittance from one electrode to the other one [5]. 
These studies are of utmost importance for nano objects, 
for which the properties change not only from object to 
object, but also on the sub-object length scale. We here 
extend the capabilities of STM further by imaging the 
inelastic response of electrons scattered off nano objects. 
As a model system we investigate dimers of dichloroben- 
zene adsorbed on Au(lll) at 5 K. 

The dichlorobenzene molecular structure is imaged as 



two triangles pointing away from each other correspond- 
ing to the structure indicated by the ball-and-stick mod- 
els in [lji [13 j. Inelastic electron tunneling spectroscopy 
(IETS) maps, shown in [l^i and b, represent the second 
derivative of the tunneling current with respect to the 
voltage at a certain voltage (d 2 I/dV 2 ), which contains 
mainly the inelastic signal as detailed below. We here 
investigate the distinctive wave pattern emerging around 
the double triangle. The pattern shows three to four rings 
that follow the contours of the molecule pairs. In addi- 
tion, these rings are modulated in intensity. The wave 
patterns are reminiscent of oscillations due to inelastic 
scattering at vibrational impurities suggested by theory 

The vibrational origin of these patterns is revealed in 
the IETS spectra on a single molecule ([2| . On the surface 
(spectra 1 and 8) only the characteristic phonon peak of 
the surface with its typical inverted symmetry (maximum 
at negative and minimum at positive voltage) is observed 
[7]. Above the molecule (spectra 5, 6, and 7) a maximum 
is observed at 10 meV. The point symmetry of this peak 
with respect to the Fermi energy identifies it as of vi- 
brational origin. At this low energy, external vibrations 
of the molecule towards the surface are expected. The 
broadness of the peak indicates that several modes are 
involved here. A peak in the same energy range, though 
with lower intensity, is observed within the region of the 
halo. The standing wave pattern in the IETS maps is 
in fact most clearly observed in the energy region of this 
peak. First principles density function theory calcula- 
tions performed on the molecular structure reveal that 
there are rotation and tilting modes around 3 to 7 meV 
and C-Cl torsion and stretch modes between 21 and 52 
meV.[T5] All those modes have a non- negligible coupling 
strength to the surface, and it is therefore likely that 
modes of this kind are responsible for the patterns in our 
IETS maps. 
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FIG. 1: a, Leftmost panel: topography with stick- and- ball models indicating the orientation of the molecules; the CI ligands 
form the two sharper edges of this triangle. Tunneling current 60 pA at 95 mV. Other panels: IETS maps recorded in constant 
current mode; note that the intensity scale is scaled differently on the different images in order to adjust to the different 
intensities of the patterns. Tunneling current 60 pA, temperature 5 K, map size: 8.7x8.7 nm 2 . b, Zoom-in on maps at opposite 
polarity ±5 mV, 60 pA. c, Schematics of the measurement principle. 



We conjecture that the interference pattern is of inelas- 
tic origin. We explain it as arising from the superposition 
of electron wave-functions before and after scattering at 
the nano-object as sketched in Fig. In contrast to 




o 

V(mV) 



FIG. 2: IETS spectra of dichlorobenze. IETS spectra 
a recorded where indicated on the STM image b, 55 mV, 14 
pA: on the surface (1,8) close to a molecule (2,3,4) and on the 
molecule (5,6,7) with modulation of 8 mV, 361.1 Hz; vertical 
lines at ±10 mV to guide the eye. The spectra are off-set for 
clarity. 



conventional Friedel-like oscillations, however, for which 
the electrons are elastically scattered and thereby the 
incoming and outgoing electrons have the same energy, 
here scattering is inelastic. Therefore, the energy of the 
outgoing wave is reduced (or increased) by one quantum 
of energy corresponding to the energy of a vibrational 
mode. 

Interpreting the standing waves in terms of inelastic 
scattering is theoretically corroborated by calculating a 
d 2 I/dV 2 map equivalent based on the inelastic scattering 
assumption in a real space Green function approach. The 
experimental map visualizes, the inelastic signal mea- 
sured in the second voltage derivative of the tunneling 
current, d 2 /(r, V)/dV 2 , at the position r and voltage 
V . Theoretically, we relate the inelastic signal to the 
energy derivative of the local density of electron states, 
dN(r,Lj)/duj. The local density of electron states (DOS), 
7V(r,c<j), represents the quantum mechanical probability 
for adding or removing an electron with energy w at a 
point r in real space. More details about the compu- 
tational method are given in the Supplementary Infor- 
mation. Calculating the local density of states and its 
derivative plays the role as the theoretical counterpart to 
the scanning tunneling microscopy technique. 

In short, we model the inelastic signal in the measure- 
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FIG. 3: Conceptual comparison of the charge den- 
sities, a, Energy dependence of the charge densities with 
and without inelastic contribution, N(t,lj) and iVo(r,Gi;) (up- 
per panel) , and their corresponding derivatives (lower panel) . 
A^o(r,cj) and its derivative have been shifted downwards for 
clarity, b, Spatial dependence of charge densities with (upper 
panel) and without (lower panel) inelastic contribution. The 
star marks the point at which the spectra plotted in panel 
a were calculated, c, d, IETS maps, i.e. dN(r,uj)/duj and 
dNo(r, uj) / duo for different energies. 



ments in terms of the real space Green function by means 
of coupling of the surface electrons to the local vibra- 
tional scattering centers in the spirit of Holstein coupling 
[T6] . The total Green function, G(r,r / ;cj), provides the 
conditional probability for removing an electron with en- 
ergy uj at the position r given that an electron with the 
same energy was created at r'. Scattering of the elec- 
tron at a local vibrational mode at the position R by the 
electron is picked up by the Green function self-energy 
E(R, cj), which provides satellite signatures in a narrow 
energy range around the inelastic modes ±co>o centered 
symmetrically around the Fermi level Ep, corresponding 
to the experimental observation. 

The central concept of the IETS measurement, 
<i 2 /(r, V)/dV 2 or dN(r,uj) / duo, is based on the fact that 
contributions to the tunneling conductance (DOS), which 
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FIG. 4: Comparison of calculated and experimental 
maps, a, Experimental d 2 I/dV 2 map at (5 =b 5) mV. b, 
Theoretical IETS map at 6 mV. c Line scans along c radial 
directions with respect to the molecular structure as indicated 
in a, b, from theoretical (faint line) and experimental (bold 
line) maps. Line scans are off-set for clarity. 



are slowly varying with energy will be suppressed upon 
differentiation, i.e. in the IETS signal. On the other 
hand, features which are rapidly varying with energy 
will be enhanced by the differentiation. This is demon- 
strated in [3^i, in which we plot the calculated DOS and 
corresponding IETS signal for the dressed (including in- 
elastic scattering) and bare (without inelastic scattering) 
system. Hence, although elastic Friedel oscillations may 
be recorded in the tunneling conductance (see [3}d, lower 
panel) such oscillations are suppressed in the IETS sig- 
nal since the effects related to this type of scattering are 
slowly varying with energy (see [3]i). Likewise, oscilla- 
tions recorded in the IETS signal have to be related to 
contributions in the DOS which vary on a small energy 
scale as illustrated by the sequence of IETS maps in [3]3. 
In particular, the inelastic scattering occurring due to a 
single molecular vibration provides a finite contribution 
to the DOS only in a narrow energy range around the 
inelastic mode, at which we expect the appearance of in- 
elastic Friedel oscillations. However, due to limitation of 
the experimental resolution, the recorded maps display 
an integrated average over an energy range of 10 meV, a 
range in which often more than one vibrational mode is 
found. This means that we are presently unable to pin 
down the energy of each specific inelastic mode. Despite 
this averaging, it can be observed in [T^i, that the IETS 
maps changes with energy, which is in agreement with 
the theoretical arguments. Theory thus proves that the 
observed pattern is of inelastic origin. 

We now compare an experimental IETS image to one 
calculated for the specific geometry of two dichloroben- 
zene molecules (Fig. [4]). On the gross scale, both images 
show a pattern with an envolope of a truncated rhombus. 
On a scale closer to the molecules, there is a finer pattern 
of dips and peaks. When comparing line scans taken in 
radial directions, see [4^, we find a fair agreement close 
to the molecular structure. The deviations at larger dis- 
tances are explained by mixing with waves emerging from 
the surroundings. The general trend of a regular modu- 
lation with decreasing amplitude is nicely reproduced for 
up to four maxima. The real space decay of the waves 
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emanating from the molecular structure is captured by 
[^(/colr-Rl)] 2 ~cos(2fc |r-R|)/(2fc |r-R|) for large 
fco|r — R|, where Hq~\x) is the Hankel function whereas 
fco is the wave vector corresponding to the inelastic mode 

UJ . 

Combination of theory and experiment thus resolve 
that the inelastic Friedel-like pattern originates from the 
electrons sampling the molecular excitation spectrum. 
The signature of the molecular vibration is carried over 
to the electronic density of the metal and leads to stand- 
ing waves around the scattering centers. These inelastic 
Friedel-like oscillations are spatially mapped by recording 
the inelastic signal, d 2 I/dV 2 . The sensitivity of the STM 
tip towards this signal is enhanced by functionalization 
via picking up a dichlorobenzene molecule (cf. [T71 fl8]). 

In this paper we have demonstrated a system which ful- 
fill the requirements for observation of inelastic Friedel- 
like oscillations emerging from a vibrational adsorbate 
on metal surface. The inelastic Friedel-like oscillations 
are not restricted to the specific system investigated 
here, which was also discussed in Ref. [14]. First, they 
should be observable on all surfaces that support an oc- 
cupied surface state, i.e. on all surface, on which elas- 
tic Friedel-like oscillations have been reported, see e.g. 
Refs. JT21 [19]. Second, all nano-objects that can be 
excited vibrationally should be able to scatter the elec- 
trons inelastically. This comprises certainly a plethora 
of molecules, but also metallic nano-clusters, that sup- 
port phonon modes. We propose that in future it will 
be possible to investigate the inelastic properties of such 
nanoobjects can be investigated by the method lined out 
in this article. 
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